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Abstract 

A modular application of the integration by fractional expansion method for evaluating Feynman diagrams 
is extended to diagrams that contain loop triangle subdiagrams in their geometry. The technique is based 
in the replacement of this module or subdiagram by its corresponding multiregion expansion (MRE), which 
in turn is obtained from Schwinger's parametric representation of the diagram. The result is a topological 
reduction, transforming the triangular loop into an equivalent vertex, which simplifies the search for the MRE 
of the complete diagram. This procedure has important advantages with respect to considering the parametric 
representation of the whole diagram: the obtained MRE is reduced and the resulting hypergeometric series 
tend to have smaller multiplicity. 
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1 Introduction 



The integration by fractional expansion technique (IBFE) constitutes a rather simple method that allows 
to evaluate Feynman diagrams starting from the corresponding Schwinger's parametric representation, and 
whose main advantage is to transfer the complexity of the direct calculation of multidimensional integrals 
to the evaluation of linear systems of equations that the method generates. This technique is particularly 
useful when the propagator exponents are arbitrary, and the general solutions obtained are sums of generalized 
hypergeometric functions, whose number of variables depends on the number of invariants in the diagram and 
also on the topological family to which the diagram belongs. 

As is the case with any integration technique, although IBFE can be applied to any diagram, it is not always 
possible to obtain simple solutions due to the number and multivariability of the resulting hypergeometric series. 
Nevertheless, for certain families of diagrams its application is particularly simple and optimal. In a previous 
work [2] we showed that a modular application of IBFE is a tool that simplifies considerably the complexity 
of the solutions : specifically we considered bubble-type modules or subgraphs. Here we extend this modular 
application to cases in which the modules are triangular loops, which then provides a rather simple procedure 
to find the solutions to a large number of Feynman diagrams. 

This work is divided in the following way. In Sec. // we apply the IBFE technique to the massless triangular 
module, and find its representing multiregion expansion for the case of off-shell external lines. Here we define 
the one-loop function Gt , which allows to describe the multiregion expansion (MRE) of this module as the sum 
over this loop function and a vertex, thus effectively reducing the loop to a node. One of the main characteristic 
features of the Feynman diagram general solutions, when the propagators have arbitrary powers, is that they 
are given by sums of multivariable hypergeometric functions. The arguments or variables in these functions 
are ratios between invariants associated to the graph, and given the structure of the resulting hypergeometric 
function, the convergence condition requires that at least this ratio of invariants be less than unity. Therefore 
from a particular diagram MRE different solutions can be extracted, each of them associated to a particular 
kinematical regime, where the differentiating element is determined by the ratio of two invariants. This form of 
differentiating the solutions is direct for certain families of diagrams [2j. Once all terms associated to a specific 
region of interest are found, summing them will give us the desired solution. Nevertheless, for topologies that do 
not belong to the previously mentioned family, and which is the most general case, this criteria for finding the 
correct terms of the solutions extracted from the MRE actually fails, because in the general case the application 
of IBFE can generate solutions which contain hypergeometric functions with a particular argument, the unity. 
Given the structure of the diagram solutions, this value corresponds to the ratio of two invariants (A/B) or 
(B/A) , when A = B. Nevertheless, both ratios are associated to different regions (different solutions) and the 
fact that they are equal does not allow to differentiate them. Here we present a method that solves this problem, 
and apply it to the specific case of diagrams that contain triangle subdiagrams. 

In Sec. /// we evaluate two different two-loop diagrams which contain the triangle subgraph in their topology. 
The first corresponds to a correction to a three-point function, with six propagators. We find the general solution 
for arbitrary propagator powers, and then particularize it to unity. In order to validate this technique, we 
compare the solution with the one obtained from a conventional method IBP (integration by parts) [4j. The 
second diagram is a self-energy correction with five propagators, for which we present a new solution for the case 
or arbitrary values of the propagator powers. Finally in Sec. IV we describe some cases in which this massless 
triangle function can be used and the generalization to massive loops. 

2 Triangular module 

In this section we will find the MRE associated to the one-loop triangular subgraph or module, considering 
a case in which the propagators are massless and the external lines are off-shell. From the MRE of this graph we 
will define the one-loop function Gt, which will allow to represent this module as a particular sum of three line 
vertices. This formulation will then be applied as examples to a self-energy and a three-point function diagram, 
both with two loops. The self-energy example will show quite clearly the advantages of this modular procedure. 
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2.1 Obtaining the MRE of the triangular module 



We start from the momentum integral representation for this graph, which is given in the massless case and in 
D dimensions by: 



G = 



(1) 



When the propagator powers {ui} are integers, this expression can be easily reduced to simpler topologies using 
the integration by parts technique [4]. The situation is more difScult for arbitrary indices {ai}, which we are 
considering in this work. Graphically the integral can be represented by the following diagram: 




(2) 



The IBFE integration technique is applied to Schwinger's parametric representation of the diagram, which can 
be obtained from the matrix of parameters [3] associated to Eq. ([1]). In the present case this matrix can be 
easily obtained and it is given by: 



Xi + X2 + X3 Xi + X2 X2 
M = I Xi + X2 Xi + X2 X2 

X2 X2 X2 

which allows to express Schwinger's parametric representation as: 

CiipI + 2Ci2 Pi.p2 + C22 pI 



(3) 



r(ai)r(a2)r(a3) 



u 



JjD/2 



(4) 



where dx — 



ai— 1 02 — 1 03 — ! 



dxidx2dx3 and U = {xi + X2 + X3). The coefficients Cij can be obtained from 



the determinants of the submatrices of the matrix of parameters (jS]): 
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C22 — 



Xi + 


X2 + 


X3 


Xl ^ 


- X2 


Xi 


+ X2 




Xl ^ 


- X2 


Xl + 


X2 + 


X3 


X2 




Xl 


+ X2 




X2 




Xl + 


X2 + 


X3 


X2 






X2 




X2 





X3{xi + X2), 



X2X3-, 



(5) 



X2(xi + 2:3). 



Replacing in (g]) and remembering that p\ — (^1+^2)^ = Pi + 2pi.p2+P2' 8'^^ after a little algebra Schwinger's 
parametric representation of ([T]): 



G 



n-=ir(a,) 



d X 



( X3X1 2\ / X1X2 2\ I 

exp y —Vi ) exp y —vi ) exp y 



X2X3 2 



(6) 



Now we will deduce the MRE of the parametric integral For this purpose we expand the exponentials that 
are present in the integrand, which leads us to the following multiple series: 
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n;=ir(a,)„^^ 



y: 0».,..,«3 {pt {pir {pir / 



"1 +"2 „«2 +"3 ^"1 +"3 



jy'-y+ni+n2+n3 



(7) 



and this allows us to obtain the MRE of the polynomial U — [xi ~\- X2 +2:3). For the denominator in the integral 
we have that: 



1 



^ ^ 4'ni,..,ne -^1 ■^2 ■^3 
ni,..,ne 



(y + iT-i + n2 + + 714 + "-5 + "e) 

r(f +771 +712 + 713) 



(8) 



and then, integrals already separated in factors of the form J dx x°'^'"^^, can be replaced by the equivalent 
constraints {a + ...). This is the final step of the fractional expansion process, and we have finally obtained the 
MRE of the diagram from its parametric representation 



G 



(-1) 



-D/2 



11^ = 1 r(aj) ni,..,n 



{pir ipir (pir 



n 



r(-J + Til + 772 + ns) 



(9) 



Where we have summarized the constraints {Aj} in the following identities: 

Al = (-J + 77l + 772 + 773 + 774 + 775 + "e) , 
A2 = (fll + 71l +772+ 774) , 
A3 = (a2 + 772 + ns + 715) , 

A4 = (03 + 771 + ns + ne) . 

From expression ^ we see that the composition of sums and Kronecker deltas (« Afc) of the MRE of this 
module or graph is 6S/4(5. 



(10) 



2.2 Definition of the one-loop function Gt 



In order to generate a systematic reduction procedure of interior triangle subgraphs of a more complex topology, 
we define from Eq. ([9]) a loop function denoted Gt, such that ([9]) can be written as: 



G= ^ GT(ai, 02, 03; 771, 772, 773) 



1 



1 



1 



ipi)-' (Pi)-^ ipi)-' ' 

where the one-loop function Gridi, 02, 03; 771, 772, 773) by comparison is given by the expression: 



GT(0l,a2,a3;77l,772,773) = — 3 — — — - J2 

l[i=i J- m) "4--" 



n 



r( -J + 771 + 772 +773) 



(11) 



(12) 



Let us look at the propagator structure of the Eq. (fTT|) . It is clear that it can be represented pictorially as a sum 
whose argument contains the original diagram with the loop reduced to an effective vertex where the external 
lines reach, as inverse propagators. The cost of this topological minimization is a multiple sum. Graphically the 
representation of the MRE of the triangular module is given by: 
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Therefore we have shown that it is possible to reduce topologies which contain this type of vertices, just as was 
done with the reduction of propagator bubbles [1] . 



2.3 Method for differentiating series of unitary argument 

Before going into the actual application of IBFE to diagrams with triangle subgraphs, it is necessary to implement 
a method that allows to discriminate whether a unitary argument in the series that comes from the MRE, 
corresponds to the limit of the ratio of invariants in a particular region or to its complementary region. This 
is based on multiplying the external momenta by certain constants (fictitious invariants) Ai [i = 1,2,3). Once 
the solution to the diagram is found, the arguments of the obtained hypergeometric series correspond now to 
ratios of these constants, allowing to separate the solutions associated to different regions of the hypergeometric 

functions whose arguments are ( -7— ) and ( — p ), with i ^ k. Once these structures have become separated. 



we take Ai = A2 = A3 = 1, and the correct solution is obtained. This simple method solves the problem 
of summing hypergeometric series, with the same argument, that actually belong to different solutions of the 
original problem. In the triangle case, we transform the external momenta in ([9|) as follows: 

pj =^ ApI (14) 
This transformation allows to obtain the modified expression of ((T3)) : 




E {A.riA^riA,) 



GT(ai, 02, 03; ni, 712,713) 




3 Applications 

The topological reduction formula (fT5|) . which corresponds to an MRE of the triangular module, will be 
useful for finding the MRE of generic graphs that contain in their geometry one or more of these triangular 
modules. We will now explain how to use formula (|15p in specific cases, showing the advantages of this procedure. 

3.1 Example I : Two-loop triangle 

As the first example of the modular reduction technique and of the topological formula (|15p . let us consider a 
diagram g with three external lines, two loops and six propagators. We take a massless theory, with mass shell 
external lines Ki and K2 {Kf = iff = 0). The diagram is: 
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3.1.1 Getting the diagram MRE 

The first step in order to find tlie diagram MRE witli tlie modular procedure is to reduce the subgraph associated 
to the indices {01,02,03}, using formula (jlSp . In this way we obtain an expression where we have eliminated 
one of the loops of the diagram, which has been replaced by an effective vertex where the three propagators 
associated to the indices {— ni, — ri2, —^3} meet. In detail: 
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Y: {A,r iA2r i^sr Criai^a^, 03; rii, ^2, ^3) 




(17) 



In Eq. P?)) we can see that one of the external lines, the one associated with the momentum Ki, is affected by 
the modular reduction on g. It is possible to actually extract this propagator => , ^^s^-„-^ from the diagram, 

since it is not involved in the loop integration, and write it as a factor that multiplies this graph. Looking at Eq. 

and its pictorial equivalent (|13p . we conclude that it is possible to eliminate from the graph the information 
related to the index of summation ni. Rewriting (|17p we thus have: 

3= E (v4i)"M^2r(A3r(x2)"^ GT(ai, 02,03; ni,n2,n3) 




(18) 



The resulting triangle allows us to apply the reduction formula (fT3|) . We then get the MRE of the topology g 
as: 
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E (Ai)"^ {A2r (^3)"^' {KlY' X GT(ai,a2,a3;ni,n2,n3) 

ni,..,n3 



X GT(a6 - ^2, as, 04 - rig; 717, ng, rig) x 




(19) 



Since each function Gt is associated to six summation indices, in order to use the same notation for these indices 
{n}, the second function Gt starts with the summation index ri'j. Therefore we get: 



5= E {A,r {A,r {A^r^ {KiY'^'^' {Kir {KiY' 

ni ,..,713 

GT(ai,a2,a3;ni,n2,»T-3) x Gria^ - 712,05,04 - n3;n7,ns,ng), 
where after replacing the functions Gt by their respective MRE, we obtain the expression: 



n^.r A, 



+ + '^2 + W3)r(-f + ny + ns + ng)r(a6 - n2)r(a4 - na) 



We have defined the factor: 



F 



(-1)" 



r(ai)r(a2)r(a3)r(a5) 



The constraints {Aj} are given by the following identities: 



(20) 



(21) 



(22) 



(23) 



Ai = (-J + ni + n2 + ns + 714 + n5 + , A5 = (-j + + rig + + riio + nn + 7212) , 
A2 = (fli + ni + 712 + "-4) , Ag = (04 - 713 + 717 + ng + riio) , 

A3 = (02 + 712 + "3 + "•s) , A7 = (05 + 77-8 + Tig + Tin) , 

A4 = (03 + Til + 713 + 716) , Ag = (ae - 712 + 117 + ng + '^-12) ■ 

FinaUy, applying the conditions of the problem {K\ — — O) , the summation indices ni , 717 and ng must 
vanish in order to have a nonzero solution. Replacing these values and eliminating the sums in (j2ip . we have 
obtain the MRE of the diagram. 

After renaming the indexes {71} =^ {1} so that they become consecutive, the multiregion expansion is finally 
given by: 



5= F Y.^lu..M (^2)'M^3)'M^3 

{1} 

n ■=! A, 



2\'6 



(24) 



r(f + ^1 + /2)r(f + h)T{aQ ~ h)T(ai ~ h) ' 



and now the constraints are given by: 

Al = (f + /l +/2 + ^3 + ^4 + ^5) , 

A2 = (ai +/i + /3) , 
A3 = (a2 + ^1 + /2 + ^4) , 

A4 = (03 + ^2 + ^5) , 



A5 = (f + ^6 + ^7 + ^8 + ^9) 

Ae = (a4 - ^2 + ^6 + ^7) , 
A7 = (as + ^6 + ^s) , 
Ae = (ae - + h) ■ 



(25) 
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Therefore the muhiphcity of the hypergeometric series that can be extracted from the MRE is one hyper- 
geometric series of the type and at most C| = 9 series of this type whose argument for this case is 
one. 

3.1.2 General IBFE solution 



Now we proceed to find the general solutions associated with the diagram coming from the MRE which 
represents it. We maintain all the indexes {aj} with arbitrary values, and for notational simplicity make A2 = A 
and A3 = B. 

The solutions we present here correspond to two regions that can be differentiated according to the ratio of 
the fictitious invariants used in order to separate them: the region where A > B and the region where B > A. 
Both are solutions of g, related by analytical continuation. In particular, in the limit A = B = I they become 
identical, as will be seen later on. 



Analytical solution in the region 

indexes {aj}: 



< 1 Let us first define the following simplifying notation for the 



- ak ■ 



(26) 



Moreover, since in practice each constraint (...) eliminates a sum in the MRE, the remaining or free sum at 
the end of the elimination process generates the corresponding hypergeometric representation qFq-i which we 
will identify as Gj , since it is the contribution obtained when the summation index rij is free in the multiregion 
representation of diagram g. 

With these definitions we can write the solutions for this region as: 



g^g (I) = Gi + G5, 
where the terms inside the sum correspond to the following functions: 



(27) 



Gi 



0-0123456 r(ai23 - y)-^("i23456 - D)^[^ - ai'i)^[^ " ai2)r(Y - a56)r(£) - 012345) 
V{a2)T{a3)T{a5)T{D - ai23)r(ai234 - f )r(^ - ai23456) 



and: 



ai23 



D 
2 ' 

_ D 
2 



D 
2 



^56 



ai2 - ^ 



(28) 



G5= [-1)-^{K?- 



Z5-ai23456 r(ai2 - ~ "2)r(Y - ai3)r(D - ai2346)r(ai23456 ^ D)T{D - Q1256) 

r(ai)r(a2)r(a5)r(Z? - ai23)r(a34)r(2^ - 0123456) 



-D — ai256, 0,3, "T ~ 0^2 



034, ^ + 1 - a 
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(29) 



Analytical solution in the region 

is given by: 



< 1 Analogously we have that in this fictitious region the solution 



9 = 9(i)=G2 + G3 + Gr, 



(30) 



with: 
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D 1 ^-«12 



G,= {-l)-^{Kl) 



r(ai23 - f )r(f - a23)r(f - ai3)r(ai23456 - D)T{D - ai2356)r(£' - ai2346) 

r(ai)r(a2)r(a5)r(a4)r(I? - ai23)r(^ - ai23456) 



J- — -T + «23) ~ T + "12356 



(31) 



^ ^ i\-r> /r^2\^-«i23456 r(a23 - t)^(t ^ «13)r(— - a2)r(ai23456 - -D)r(— - ai5g)r(D - ai2346) 

G3= (-1) (Asj — - — — FiTFTsS : ^ 



r(a2)r(a3)r(a5)r(I) - ai23)r(a234 - D)T{^ - ai23456) 



' T + l-a23, l+ai56-T 



and finally: 



(32) 



Gr= (-l)-^(i^l) 



0-ai23456 r(a23 — ■j)r(-^ - ai3)r(-^ — a2)r(ai23456 - -D)r(-j - ai56)r(_D - ai2346) 

r(a2)r(a3)r(a5)r(D - ai23)r(a234 - D)T{^ - ai23456) 



+ 1 — ai23456i Y ^ '^56' ^ ^ '^1256 
D + 1 — 012356, T + 1 ~ ^^156 



(33) 



3.1.3 Particular case : unit indexes 

In the actual Feynman diagram evaluation, in general the corresponding momentum integral has indexes or 
propagator powers that are equal to unity. From the previous results, making = 1 (j = 1, 6), we obtain 
the solution: 



Analytic solution in the region 



where now: 



< 1 The solutions in this region are given by: 



(34) 



Gi = (-1)-^ {Kl) 



D-6 r(3 - f )r(6 - D)r{^ - 2fr{D - 5) ^ 
r(D-3)r(4-f)r(3i^-6) 



(35) 



and: 



Gs = (-1)-^ {Kif-' ^C^-TnT - l)r(f - mo - 5)r(6 - D)T{D - 4) 



D-A, 1 
2 



(36) 
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< 1 In this region the solution is: 



Analytical solution in the region 

5 = 5(1) =G2 + G3 + G7, 
where the hypergeomctric representations 62,3,7 are now given by: 

n_6 r(3 - f )r(f - 2)2r(6 - D)riD - s)^ 



G, = (-1)-^ (/^r 



r(z?-3)r(^-6) 



A^-' 2F1 



0, 1 

6- D 



or equivalently: 



D-e r(3 - f )r(f - 2)^r(6 - D)T{D - 5y 
r{D - 3)r(^ - 6) 



X A~ 



(37) 



(38) 



(39) 



and finally: 



Ga = (-1)-^ {Klf-' ^C^-^n Y 2)r(f - l)r(6 - D)TiD 5)r(f - 3) 



1 ^ - 2 



r(i? - 3)r(^ - 6)r(3 - f ) 

1 



D- 



-e r(2 - f )r(f - 2)r(f - i)r(6 - z?)r(f - 3)r(i? - 5) 



T{D-3)Tm-6) 



D~5 



(40) 



(41) 



3.1.4 Analytic continuation and dependence between the solutions obtained with IBFE 

It is possible to show that in the limit A = B = 1 the solutions ([34|) and ((37|) are equal, and therefore both are 
related by analytic continuation. For this purpose we use the following identity: 



2F1 



a, b 
c 



1 = 



T{c)r{c - a-b) 



r(c-a)r(c-6)' 

Using this formula we easily find that the summed series have the following form: 

-Drr.2^D-e r(f - 2)^r{D - 5)r(5 - D) 



r(D-3)r(^ -6) 



G5 = (-l)-^ {Kl) 
On the other hand we have that 



^ , D-6 r(2 - f )r(f - 2)r(f - i)r(D - 4)r(i? - 5)r(5 - d) 



r(D-3)r(^ -6) 



G, = (-1)-^ {Klf-' ^^^^^ ~ ^^'^^^ " ^^^^^ " 



and : 



r(D-3)r(^ -6) 



G ( 1)-^ (K^)''-' ^^^^^ ^^^^^ ^^^^^ ~ ^^^^^ ~ (4 - f ) r (5 ^ D) 

^ ^ > ^ 3j r(i?-3)r(^ -6)r(3- f )2 



(42) 



(43) 



(44) 



(45) 



(46) 
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The hypergeometric function contained in the term G-j and which has the form i_Fo: can be transformed in 
another one of the type 2^1, and then we can apply the identity ([1^ . 



and then we get that 



iFi 



G7 = 0. 



a, P 



(47) 



(48) 



After a bit of algebra, and using known Gamma function properties, it is straightforward to show that the 
solutions ([M)) and (P7|) are equal in the limit A = B = 1, giving 

(Gi+G5)-(G2+G3) = 0. (49) 
Therefore, defining Gibfe as the solution for the graph in (|16p . the final solution in general can be written as: 

G2 + G3 



G 



IBFE 



= < 



(50) 



Gi + G5. 



3.1.5 Validating the solution Gibfe- Comparison with a conventional method : Integration by 
parts (IBP) 

When the propagator powers are integers, it is possible to evaluate such topologies with the integration by parts 
method |4]. In the particular case of diagrams that contain subgraphs, this technique can be quite useful in 
order to simplify the integration problem, since it allows to simplify the geometry of the diagram. 

In the particular case of the diagram which contain triangle subgraphs, this technique can be in order to 
simplify the integration problem. One very useful formula, which can be deduced applying IBP to a triangle 
graph, and which is known as the triangle identity, is the following: 




1 



(D - ai - a2 - 2a3 



\ 




(51) 



+ 02 



where it has been assumed that the indexes ai,a2, and 03 have arbitrary values. We now apply directly this 
identity to the diagram depicted in (fTH]) . but taking all the propagator powers to be unity. The result that we 
obtain is: 
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The third diagram in the left-hand side of the equation vanishes (external on-shell line condition). The second 
diagram can be further simplified using again (|5ip . giving: 




G(ai, 02) 



1 



in 



(54) 



a2 



where: 



G(ai,a2) = (-1) 



-D/2 



r(ai 



a2 



#)r(f -ai)r(f 



02) 



r(ai)r(a2)r(i^-ai -02) 



(55) 



In the case of the triangle with conditions Pi — P2 — one has the fundamental formula: 





P3 = G (01,02, as) {P3) 



2\ -3-— 11— 13— 02 



(56) 



where the factor G(ai, 02, 03) is given by: 
G(ai,a2,a3) = (-1)"''/' 



r(ai + 02 + 03 - -j)r(-^ - a2 - a3)r(-| - oi - 02) 
r{ai)T{a3)r{D - oi - 03 - 02) 



(57) 



Applying formulae ([M)) and ([55)1 , the solution for the diagram is then: 
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g = G 



IBP 



D-4: 
1 



2G(1,2)G(4 - f ,1,1) - -^G(2,2)G(1,4 - f ,1) 



(58) 



D-5 



2G(1,3) [G(l,4-f,l)-G(l,l,4-f)] 



Choosing one of the solutions for Gibfe ([50| . it can be easily shown that the solution of this diagram evaluated 
using the IBFE technique gives the same result as with IBP: 



G 



IBFE — GjBP 



(59) 



This is then a concrete confirmation of the IBFE technique, applied modularly to a diagram with triangle 
subgraphs. The important point is that this has allowed to obtain the MRE of the graph in a noticeably more 
direct and systematic way than finding it by taking Schwinger 's parametric representation of the whole diagram. 



3.2 Example II : Two- loop propagator correction 

A more relevant example is the two-loop propagator correction, with five internal lines. 



G = (x> 




(60) 



3.2.1 Obtaining the diagram MRE 

In order to find the diagram MRE we apply formula (jl5p on the left-hand side of the diagram, obtaining the 
following graphical expression: 



G= J2 (^i)"M^2)"MA3)"^GT(ai, 02, a3;ni, 712,^3) 




Using then the massless bubble one-loop functions Ga |I] , we quickly obtain the diagram MRE as 



G 



{Ax) '(^2) M^s) GT(ai,a2, as; ni,n2,n3) X Ga (04 - 711,05 - »^3;»^7) 

ni ,..,r!3, ii 

(-^2 - n-j) 
® * ® 



or equivalently: 



(61) 



(62) 



G = ^ (^1)"' (^2)"' (^3)"^ GT(ai,a2, 03; 711,712, "3) x GA(a4 - 711,05 - 773; 717) ^ 



The one-loop functions Gt and Ga a-re given by the expressions: 
GT(ai, a2, as; 77i, 772, 773) = — — — ^ ^ 



(p2 



n =1 A, 



n=ir(a,) 



1 (y + 77l + 772 + 773) 



(63) 



(64) 
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whose constraints are: 




(65) 



and for Ga we have: 

Ga (04 - "1, as - nsinr) = 
where the constraints become: 



{-l)-D/2 



r(a4 - ni)r(a5 - n^) 



E 



j=t> "J 

r(f + 717) 



(66) 



A5 = (-f + + ng + rtg) , 
Ag = (04 - ni + 717 + ns) , 
A7 = (as - 713 + 717 + ng) . 

With this information we finaUy are able to write the diagram MRE as: 



(67) 



G 



rii^l r(aj) ni,..,ng 



{A,r iA2r w {p') 



"2+17 



r(f + 711 + 712 + 7i3)r(f + 717) r(a4 - 7^l)^(a5 - 7^3) 



(68) 



This expansion contains OS and 76, and therefore the solution of this diagram corresponds to a double series and 
since we are dealing with a propagator, the kincmatical variable (p^) that is present does not appear explicitly 
as an argument in this series. We expect that for this type of diagrams the solution will not be a single term, 
since it does not correspond to an optimal topology (in the sense that the multiplicity of the resulting series 
does not depend on the number of loops or equivalently only depends on the invariants of the physical process 
i). 

In order to simplify the notation in ((68|) we take Ai = A, A2 — B and ^3 = C. With these invariants six 
regions of interest can be generated, each one with a different solution, but in the limit A = B = C = I they 
are all identical. These regions are given by: 

A<B <C, 
A<C <B, 
B <A<C, 
B <C <A, 
C <B <A, 
C <A<B. 



(69) 



We arbitrarily take the region in which {A < C < B). Then when we extract the different terms of the MRE 
(bivalued hypergeometric functions) we only sum those that have as one of their arguments the following simple 
combinations: 



or equivalently all the double combinations generated when combining these simple combinations: 
(4' ■§) ' Given these conditions, the general solution that we get starting from is the following: 



G ^ gi + 92 + 93 + 94 + 95 + 96 



(70) 



where we have that: 



14 



( 2xc-«i2345 r (ai3 - f ) r (f ~ a23) r (f - ai) r (01345 - g) r (f - as) r (ZJ - a^) 

^ ^ r(ai)r(a3)r(a5)r(i?-ai23)r(ai34-f)r(^-ai345) 



2:2:2 
X i^2:l:l 



{a2, §-ai}, {1- as, -J -05} 

{^-ai345, £' + 1-01345} , {l-f+a23}, 

{§ + l-ai34, £'-0134} £ ^ \ 
{f + l-ai3} B' B J' 



(_-,\-D , 2\-D-ai2345 r ("J ~ '^^s) r (£) - 02345) r (012345 — -D) r (-j - 013) r (0345 — ^) 
^~ ^ ^ r(ai)r(a2)r(a3)r(f)r(i?-ai23) 



C^2-l-l ( «2345, § - 0145} , {1 - 04, § - 04} , {ai2345 " D, 0345 - 

■■ I {|_a4, 1-04}, {l + ai3-f}, {f} 



A C 
C'B 



^_^yD ^^2^^-^12345 r (ai3 - -j) r (ai2345 — -D) F (D - ai34s) r — 2a3 — 01245) r (0345 



r (ai) r (02) r (03) r ) r (z? - am) 



p2:l:l {-O - 0134S, ^ - 2a3 - 01245} , {f + ^ ~ '^134> - O134} , 

{£) - 0134, f + 1 - 0134} , {§ + 1-013}, 



{^12345 — D, 0345 — y} 
{f} 



A C 
C'B 



/ 2\^-«i2345 r (f - 013) r (f ^ "23) r (0123 - f ) r (04s - f) r (f - 05) r (f - 04) 
^ ^ r(oi)r(o2)r(os)r(i?- 045)r(£)-oi23) 



XF2-l-l( {^3, ai23-f}, {1-04, f- 04}, {1-05, f- OS 

■ ■ V {D-a45, 1 + f -045}, {l-f+oia}, {l-f + a23} 



A C 
B'B 
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^ ' ' r(ai)r(a2)r(a3)r(2D-ai245 - 2a3)r(ai34-# 



and finally: 



56 



0235) 



X 



.3r-- 



ai245 + 2a3 - ^) ^^_„^3^,3_«^^_„^3 

r ("235 - #) 



2:2:2 
X F2:l:l 



{£) - ai23, -7 - 03} 
{2i?-ai245 - 2a3, ^ + 1 - 2c 

{f + : 



{f + 1 



)j 5 12'" fl235, D — 0235} , 

2a3 — ai245 } , 



{# + 1-023}, 



- 1 - ai34, D - ai3i} 
{t + 1-«13} 



C A 



X B-°iC*- 



ai3) r (a23 - t) r (f - «2) r (t - "4) r (02345 - iJ) r (d - 0235) 
r (02) r (03) r (04) t{d- 0123) r (0235 - f) r - «2345) 



(123 



2:2:2 
X F2:l:l 



{01,-7-02}, {1-04,^-04}, 

{^ - a2345, l + D- 02345} , {1 - f + 013} , 

n _ -1 . \ 



{1 + § - O235, -D - a235} 
{l + T-«23} 



(75) 



(76) 



■[i + - - 023) B B J 

The previous results have been written in terms of the functions F and F (for more information see Ref. [2] ), 
which are defined as: 

p r u 



l-si" / {7i,--,7p} {ai,..,Or} {ci,..,Cu} 



p:r:u 



p r u 

n(7,) n — m n K)n n 



and similarly we also have: 



J = l 



j )n-~-m 



n! to! ' 



j=i i=i 



(77) 



Fq^\i [ {7i:--:7p} {ai,..,ar} {ci,..,c„} 



00 n (7,) n (flj)" n (cj)m ^ 



This last series is called Kampe de Feriet function. Finally the solution is obtained by simply taking A = B = 
C= 1. 

We have solved this diagram in a modular or loop by loop form. In order to appreciate the advantages of this 
way of applying IBFE, it is important to compare the MRE (f68|) with the MRE of the same diagram obtained 
by the application of IBFE to the parametric integral of the complete diagram. For this purpose let us consider 
first Schwinger's parametric representation of the complete diagram, which is given by: 



G 



(-1) 



-D 



exp 



n r(o,) ^0 



d X 



(xi + xi) [x'i + a:;4 + ^5) + •^3 (2:4 + 2:5) 



\{xx + 2:2) (x3 + 2:4 + xg) + a;3 (x^ + x^)\ 



(79) 
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where the polynomial F is given by: 



F = [xiX2 {x3 + X4 + X5) + X4X5 {xi + X2) + 2:30:40:5 + 0:12:30:5 + 0:22:32:4] . 
If we now obtain the diagram MRE using ([79]) we can compare {Table I) this MRE with MRE ([68| : 



(80) 



IBFE complete diagram IBFE modular 



Multiplicity multiregion series (tr) 
Kronecker deltas of the expansion (i5) 
Multiplicity resulting series (cr — 5) 
Possible contributions to the solutions (C^) 



12 9 

8 7 

4 2 

495 36 



{Table I) 



We find that the MRE obtained in terms of the one-loop functions Gt is considerably reduced with respect to 
the one that we get when IBFE is applied to both loops of the diagram simultaneously, not only in the number of 
sums and Kronecker deltas, but also the difference between them is reduced, which implies that the multiplicity 
of the resulting series is less, a nontrivial fact that simplifies considerably the complexity of the solutions. 

4 Commentaries 

The calculation technology that we have presented here can be useful in order to find minimal expressions 
for the MRE of diagrams that contain triangular subgraphs, such as those that we see in the following figure: 



The loop function Gt is very useful since it opens up the possibility of evaluating other families of graphs and 
to obtain solutions with less complexity and less number of terms when compared with the application of IBFE 
to a parametric integral of the complete diagram. Although we have presented here the one-loop function Gt 
for the massless case, it is certainly possible to consider also massive propagators and define new loop functions 
for these cases. It will be always possible to use here the same reduction procedure for diagrams that contain 
bubbles or triangles that we have used for the massless cases. 

The same idea can be extended to consider one-loop functions for modules that contain four or more prop- 
agators. Nevertheless, the reductions are no longer so simple as for the bubble and triangle subgraphs. The 
important point is that the modular IBFE application, particularly when it is done loop by loop, allows to obtain 
a minimal representation for the diagram MRE, from which the analytical solution follows. 

5 Conclusions 

The general solutions that can be obtained in the evaluation of Feynman diagrams always correspond to 
multivariable hypergeometric series, whose multiplicity can be deduced directly from the MRE that represents 
the diagram. One of the characteristics of the IBFE technique is that it has a lower bound that fixes the 
multiplicity /i of these hypergeometric series, and which is related to the minimal number of invariants {N) that 
characterizes a specific Feynman diagram, described through the formula: 




etc. 



ti>{N-l). 



(81) 
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The greater-equal sign is here included since besides the number of invariants one must consider also the topolog- 
ical family to which the diagram belongs. If this is not the optimal for IBFE application, in the sense described 
in [2], then the multiplicity is also increased according to the number of loops in the diagram. This multiplicity 
can be directly recognized from the MRE of the diagram, and it is equal to the difference between the number of 
sums and constraints (Kronecker deltas) that are present. Solutions with multiplicity (which in reality is not a 
series, and contains only one term), or multiplicity 1 or 2, correspond to series whose properties are extensively 
discussed in the literature O O [71 [5] ■ Nevertheless, for triple or higher multiplicity series, the information is 
rather limited and it is therefore important to have methods that allow to reduce the diagram MRE, not only 
in the sense of lowering the number of sums and deltas, but also to reduce their difference, which is something 
that we have shown here can be done modular ly (see Example II). 

The idea of the modular application of IBFE can be easily extended to one-loop modules or subgraphs 
that contain four or more propagators. The result will be in general a reduced MRE, that will also contain 
simpler solution, i.e. hypergeometric functions of lower multiplicity. In the same way, the idea of using fictitious 
invariants for selecting the correct solution associated to a diagram can be also generalized to all the cases in 
which hypergeometric functions with unit arguments are generated. 

The modular application of the integration technique IBFE has shown to be a powerful calculational tool, 
not only for its simplicity but also because it allows to calculate many families of graphs. Here we have seen that 
the systematical graphical procedure that we have implemented for triangular modules, makes the reduction 
just as simple as in the case of bubbles contained in a diagram. 
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